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ABSTRACT 



We propose a method for recovering the shape of the mass power spectrum 
on large scales from the transmission fluctuations of the Lyman-alpha forest, 
which takes into account directly redshift-space distortions. The procedure, in 
discretized form, involves the inversion of a triangular matrix which projects the 
>- ■ mass power spectrum in 3-D real-space to the transmission power spectrum in 

1-D redshift-space. We illustrate the method by performing a linear calcula- 
tion relating the two. A method that does not take into account redshift-space 
anisotropy tends to underestimate the steepness of the mass power spectrum, in 
the case of linear distortions. The issue of the effective bias-factor for the linear 



distortion kernel is discussed. 
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■ 1. Introduction 

. P.. 

In an elegant paper, Croft et al. ( |1998|) introduced a method for recovering the shape 
of the three-dimensional primordial mass power spectrum on large scales from the one- 
dimensional transmission power spectrum of the Lyman-alpha forest. They observed that 
the two are related by an integral of the form: 

*■<*•> «JC^ (1 > 

where fen is the wave-vector along the line of sight, k is the magnitude of the three-dimensional 
wave-vector, and P and P are the one- dimensional redshift-space transmission power spec- 
trum and the three-dimensional redshift-space mass power spectrum respectively. It was 
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suggested that redshift distortions merely change the normalization of P from its real-space 
counterpart, and so a simple differentiation of P would suffice in recovering the shape of the 
three-dimensional real-space mass power spectrum. 



Redshift distortions (see [Hamilton 1997| and references therein), however, imply that P 
is in general a function of k\\ as well as k, in which case differentiation of P alone would not 
recover the true shape of the three-dimensional real-space mass power spectrum. 

We show in §|3] how to perform the inversion from the one- dimensional redshift-space 
transmission power spectrum to the three-dimensional real-space mass power spectrum cor- 
rectly, for general, not necessarily linear, redshift distortions. It involves the inversion of a 
triangular matrix, which acts as a distortion kernel. We illustrate the method in §[| with 
a perturbative example (i.e linear distortions), and demonstrate that the method of simple 
differentiation generally outputs a real-space mass power spectrum which is flatter than the 
true one. We end with some concluding remarks in §|5|. 

Before we proceed, however, let us clarify our notation on the various power spectra 
treated in this paper. 



2. A Note on Notation 



To avoid a proliferation of superscripts and subscripts, we adopt the following convention 
for the various power spectra, P, discussed in this paper. We use ~ to distinguish between 
one-dimensional and three-dimensional power spectra: P is 1-D and P is 3-D (i.e. P has a 
dimension which is the cube-root of that of P) . To distinguish between the three-dimensional 
redshift-space (anisotropic) versus the three-dimensional real-space (isotropic) power spectra, 
we rely on either the context or explicit arguments of the power spectra: the former is denoted 
by P{k\\,k) while the latter, being isotropic, is denoted simply by P(k). In this paper, all 
one-dimensional power spectra, on the other hand, are implicitly in redshift-space. Finally, to 
tell apart the power spectrum of density from that of transmission/flux, we use superscripts: 
P p versus P* , where p denotes the density and / the transmission. 



2 Croft et al. (1998) in fact differentiated the Gaussianized transmission power spectrum rather than the 
transmission power spectrum itself. Their investigation seems to indicate that the two give very similar 
results, except that the former yields smaller error-bars. We will consider the non- Gaussianized version of 
their method in this paper for simplicity. 
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3. General Non-perturbative Formula 

The three-dimensional, generally anisotropic, power spectrum of some random field 
is related to its one-dimensional projection through the following integral ( [Kaiser fc Pea- 
cock 19911) 

00 , s kdk 



p(* n ) = J k P(h,k)— (2) 

where k\\ is the wave- vector along the line of sight, and k is the magnitude of the three- 
dimensional wave-vector i.e. k 2 = k 2 + k\ where k± is the magnitude of the wave-vector 
perpendicular to the line of sight. We assume that P is independent of the direction of k^, 
by azimuthal symmetry, as is in the case of redshift distortions. Note that we have used 
P for the three-dimensional power spectrum, to distinguish it from P, its one-dimensional 
counterpart. 

The power spectra are related to the three-dimensional, generally anisotropic, two-point 
correlation function £ by the following: 



P(Jfc||) = 2 J 0) cos(Jfe||«||) du\\ (3) 

_ /"OO /"OO 

P(k\\,k) = An / u±) cos(k\\u\\) Jo(k±u±) u±du±du\\ 

Jo Jo 

where Jo( r ) is the zeroth order Bessel function. The two-point correlation £ depends on the 
magnitude of ui but not its direction, again by azimuthal symmetry. 

The un coordinate above stands for the velocity along the line of sight (in kms" 1 ) i.e. 
un = c(A — A)/A where A is the observed wavelength, A is the mean wavelength of interest, 
and c is the speed of light. 

The u± coordinate stands for the transverse distance in velocity units i.e. u± = 
Hx±/(1 +z), where x± is the actual comoving transverse distance, z is the mean redshift of 
interest and H is the Hubble parameter at that redshift. The mean redshift and the mean 

— o 

wavelength are related by A = A a (l + z), A a = 1216 A. 



The Fourier counterparts of un and uj_ are kn and k±. Occasionally, we will abuse the 
notation by using the (un, ku) pair to denote the coordinates in wavelength units i.e. (A — A) 
and its Fourier transform. 

The effect of redshift-space distortions on the power spectrum, at both small and large 
scales, can be described by: 



P(k h k) = W(k\\/k,k)P(k) 



(4) 
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where P(k) is the isotropic power spectrum in the absence of peculiar motion, and W is 
a suitable distortion kernel. Note that we rely on explicitly displaying the arguments to 
distinguish between the isotropic and the anisotropic power spectra. 

Finally, putting eq. (|J) into eq. (pj), it can be seen that the one-dimensional redshift- 
space power spectrum is related to the isotropic three-dimensional real-space power spectrum 
by a linear integral equation: 

P(*||) = j k W(kn/k,k)P(k)^ (5) 

Thus far, we have not specified the actual random field whose power spectrum we are 
interested in. The random field could be the mass overdensity 5 = Sp/p or the transmis- 
sion/flux overdensity 5f = Sf/f, where / = e~ T , f = (f), 5f — f — f, and r is the optical 
depth. We will use P p or P p to denote the mass power spectrum and P* or P* to denote 
the transmission power spectrum. 

The one-dimensional redshift-space transmission power spectrum can also be related to 
the three-dimensional real-space mass power spectrum by an effective kernel, which we will 
call W fp : 

pf(k\\)= I W fp {k ll /k,k)P p {k)—- (6) 
-/fen Ztt 

In discretized form, this is equivalent to: 

P f = A P p (7) 

where the power spectra are represented as vectors and A is an upper (or lower) triangular 
matrix, which is invertible if none of the diagonal entries of A vanishes. The special case 
considered by Croft et al. (|1998| ) corresponds to W^ p = const., where inverting the above 
matrix equation is equivalent to the differentiation of Pf(k\\). 

The problem of eq. (0) is of course that P^ (k\\), for any given k\\, depends on an infinite 
vector: P p (k) for all fc's, from kn to, in principle, oo. To make it useful for computation, 
we have to truncate the infinite vectors somehow. Suppose one is given a finite vector of 
P^(fc||), for k^ < fey < k^ say. Eq. (|7|), in component form, can be rewritten as: 

k f 

P'(* n )-A = E A(k h k)P p (k) (8) 
fc=fcy 

where A = /°i W fp (k {l /k, k)P p (k)kdk/2-K ) and A(k\\,k) = W fp (k {l /k, k)kdk/2n. A{k\\,k) can 
be regarded as a triangular matrix in the sense that A(k», k) can be set to zero for k < k\\ 
by the virtue of the lower limit of summation in eq. (S). 
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By inverting eq. (|8]), we can in principle determine P p (k), with A left as a free parameter. 
We can do better, however, by the following observation: since P p (k) is generally a rapidly 
decreasing function of k for sufficiently high k's (~ k~ 3 , or faster if p is equated with the 
baryon density, see footnote in §4.2| ), assuming W fp (k\\/k,k) does not increase significantly 
with k, one can see that A can be made small by choosing a sufficiently high truncation k^ . 
Therefore, inverting eq. (§) by ignoring A altogether would still give accurate estimates of 
P p (k) for fc's sufficiently smaller than k^ . We will illustrate this with an explicit example of 
A or W* p in the next section. 



4. A Perturbative Example 

In this section, we will perform a linear calculation of P* , and we will assume the actual 
shape of P* on large scales, even in the presence of nonlinearities on small scales, agrees 
with that of the linear prediction, while its amplitude might not. This is in the spirit of 
Croft et al. ( 1998|) who argued that, ignoring redshift-distortions, P^ should be proportional 



to the linear P p on large scales, even though the mass fluctuations have gone nonlinear on 



small scales. The reader is referred to Scherrer & Weinberg ( 1997 ) for arguments on why 
that is reasonable, in the context of local-biasing (the mapping from 5 to the optical depth 
or transmission can be seen as some kind of local-biasing; see §||] for subtleties however). 

The output of our calculation would be a distortion kernel W^ p (eq. |J), which may 
or may not be the true kernel on large scales if the mass density field has already gone 
nonlinear on small scales. We will have some more to say about this in §[|. Nonetheless, it is 
unlikely that the true W^ p is equal to a simple constant i.e. the general method of inverting 
a triangular matrix outlined in eq. (|7]) should be used, rather than mere differentiation. The 
perturbative example set forth in this section should be seen as an illustration of the method. 



4.1. Linear Fluctuations 



To derive the linear theory limit of P* , let us start with the following general expression 
for the optical depth (see e.g. [IVliralda-Escude &; Kees 1993| ; |Hui et al. 1997] ) : 

-i 



T Mil 



E 



"HI 
l + Z 



ds 



dx 



Urydsw 



Or 



CaO" 



exp[ 



Mil 



T 



(9) 



where nm is the proper number density of neutral hydrogen, z is the mean redshift of interest, 
x\\ is the comoving spatial coordinate and the integration is done over the velocity s\\ along 



the line of sight. Velocity is related to distance by sn = x\\Hj (1 + z) + v 



pec ^|| 



where v 



pec 
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is the peculiar velocity along the line of sight, and x\\ = is the position where the redshift 
due to the Hubble expansion alone coincides exactly with z. The Jacobian \dsu/dx\\ | _1 
multiplying the proper density nni gives us the neutral hydrogen density in velocity-space, 
and the summation is over multiple streams of xy's at a given s\\. 

The thermal profile is given in the second equality, with a a o being the Lyman-alpha cross 



section constant flRybicki fc Lightman 19791) . The width of the profile is bp = \j2ksT '/m p 
where T is the temperature of the gas, ks is the Boltzmann constant and m p is the mass of 
a proton. 

Three pieces of physics remain to be specified if one were to relate the optical depth 
and the mass distribution: 1) ionization equilibrium implies that nm oc [1 + <5{,] 2 T -0 ' 7 , where 
5b is the baryon overdensity; 2) the temperature-density relation T = Tq(1 + Sb) 1 ^ 1 , where 



To is the mean temperature at 5b = and 7 is determined by reionization history ( [Hui fc 
Gnedin 1997| ); 3) the baryon distribution is smoothed on small scales with respect to the 
mass distribution (<5 b <-» 5; see below). 



Without giving further details (see |Hui fc Kutlcdgc 1997| ), in the weak perturbation 
limit (5 T <C 1 where 5 T = (t — f)/f), it can be shown that: 



W«.> = /[p-a7(7-i)]*-^ + (7-i4§£ 



W(s\\ - u\\)ds\\ 



W(s\\ 



U\\ 



>T 



exp[-(s|| -«||)7 & t ] > 



(10) 



where bx = y2kBT /m p is the thermal broadening width at temperature T , and W is 
simply a Gaussian smoothing window. 



Using the fact that 5f oc 5 T to linear order (5f = (/ — /)//, with f — e 
deduce the one-dimensional power spectrum of the transmission: 



we can 



re 

P^kw) = Aexp[-kj/kf] / 

Jkn 



[2-0.7( 7 -l)]+/nJf 



4 K \\°T 



P p {k)e 



_ k 2j k 2 kdk 



where /n = dlnD/dlna with D being the linear growth factor and a the Hubble scale 
factor (see Peebles 1980| ). The three-dimensional isotropic real-space mass power spec- 
trum is denoted by P p (k), and kp is the scale of smoothing due to baryon-pressure i.e. 
P p (k) exp[— k 2 /kp] gives the power spectrum of the baryons. As argued by Gnedin & Hui 
( |1998|) P|, kp 1 should be given by y/2H(l + z)^ 1 fj xj, where xj is commonly known as the 



3 The kp here is equal to the kp in Gnedin & Hui divided by y/2. 
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Jeans scale. The latter is equal to , ykBT /4:7ra 2 Gpfi, where /i is the mean mass per particle 
and p is the mean mass density. The numerical factor fj relating kp 1 and xj should be 0(1), 
its precise value depending somewhat on the reionization history ( |Gnedin fc Hui 1998| ), but 
it should have an insignificant effect on our work here, because we are interested primarily 
in the large scale fluctuations. 

The other smoothing scale kf, should be equal to v^/^To due to thermal broadening, 
but we can allow it to be more general to include the effect of finite resolution as well: 



Ifcff-^ b 2 ~b 2 | FWHM2 (12) 



where FWHM is the resolution full-width-half-maximum. 

The proportionality constant A for eq. (|TT|) should be equal to f 2 within the context of 
linear theory. However, in the spirit of Croft et al. ( |1998| ), we assumes the linear prediction 
gives the right shape but not necessarily the right amplitude for the power spectrum on large 
scales (see for discussions). Hence, A will be left as a free constant. 

Note that the integrand in eq. (|Tl~D is precisely of the form shown in eq. (0). In fact, 
on large scales (small k\\ as well as k i.e. small compared to fcjj, 1/&t an d k F ), modulo 
multiplicative factors, it reduces to the famous Kaiser (|1987| ) result, if one identifies [2 — 
0.7(7—1)] with the usual galaxy-bias-factor. Interestingly, the smoothing factor exp[—fc 2 //cjj 2 ] 
is exactly of the form commonly used to model nonlinear redshift distortions on small scales 
(e.g. |Peacock fc Dodds 1994[ [Heavens fc Taylor 1995| , but see also [Fisher et al. 1994| |Cole 



|et al. 19U3| ). We will take advantage of this fact, and estimate the effect of small scale 
distortions on the inversion procedure at large scales by allowing kf< to vary. 



4.2. Inversion on Large Scales 



Motivated by eq. (|H]), we consider the following inversion problem: how to estimate 
P p , on large scales, from P* , for 



pf(k\\) = / W fp (k ll /k,k)P p (k) 



kdk 
~27 



(13) 



where 



W fp (k l{ /k,k) = A'exp[-kl/kf]exp{-k 2 /k 2 F ] 

fn kj 



(14) 



7-1 



2 - 0.7(7 - 1) k 2 4[2 - 0.7(7 - 1)] 



k\b 2 TQ 
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where A' is a constant .0 

The above p is the actual distortion kernel we will use to compute Pf for some given 
input P p . However, for the inversion problem {P* — > P p ), we will not assume we know all 
the parameters in W* p , or even the precise form of W* p , except that, on large scales, it is 
equal to 

W/'(k,/k,k) = A' 



kP 2 



= 2-0.7^-1) (15) 
The j3f here is the analog of the usual j3 discussed in the context of galaxy-redshift-distortions, 
and 2 — 0.7(7 — 1) is the equivalent here of the galaxy-bias factor. Note that the scales of 
our interest are much larger than the thermal broadening width, hence the dropping of the 
term involving b\ . 

In other words, for the inversion problem, we assume we know the distortion kernel in 
the linear regime (eq. [D|, with the single parameter /3f), but otherwise do not have any 



other information regarding the full distortion kernel W^ p (eq. fl4j| ) on small scales. This 
is intended to mimic the real-life situation we find ourselves in: that we understand linear 
distortions rather well, but do not have a good grasp of nonlinear distortions on small scales. 
We will use the extra parameters in the full kernel (which we "pretend" we do not know 
in the inversion procedure) to simulate the effect of nonlinear distortions on our inversion 
procedure (in particular, the parameter kL which coincides with a factor commonly used to 
model nonlinear distortions in galaxy surveys; see e.g. |Peacock fc Dodds 1994| ). 



Let us split the integral in eq. (|13D into two parts, a part that we think we understand 
based on perturbation theory, and a part that takes care of the small scale distortions which 
we do not necessarily have a good handle on: 

rki, _, kdk r°° kdk 

pf(k\\)= Wl P (k ll /k,k)P p (k)—+ W fp (0,k)P p (k)— (16) 
Jk\\ 2n Jk* 27r 

where we have assumed: 1. W^ p ~ W{ p for k < k+; 2. fen is sufficiently small so that 
k\\/k+ ~ 0. The second term on the right then plays the role of A in eq. (^]). As explained in 
the last section, the above set-up is then suitable for an inversion analysis. One can imagine 
obtaining P p given Pf for some range of fcy's, by inverting the matrix W/ p kdk/2ii (which 
is restricted to its upper, or lower depending on one's convention, triangular entries by the 
limits of integration), treating A as a free parameter or ignoring it altogether. 



4 Note that an alternative would be to group the baryon-smoothing factor exp[— k 2 /kp] together with 
P p (k) instead of with the rest of the terms in the distortion kernel W^ p . Our inversion procedure can then 
be viewed as an attempt to recover the baryon power spectrum P p (fc)exp[— k 2 /k 2 F ] rather than the mass 
power spectrum itself P p (k). However, the two coincide on large scales. 
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Instead of doing so, we will rewrite eq. (16) into a form that is closer to the original 
analysis by Croft et al. ( |1998| ), thereby making manifest the differences from our procedure 
suggested here. 



By taking the derivative of eq. 
A'P p {k = k t ) = 



161) with respect to kn, it can be shown that 
dP f 



(1 + p f fk 



dk\\ 



k*± _, rlh 



2tt 



(17) 



k ** ~ dk 
A'pp(k)k- 3 ^ + C 2 



where we have used the form of W[ p in eq. (|T|) . The value of ki for which we will perform the 
inversion would range from some maximum k** to whatever small ki (large scale) one might 
wish. The constraint is that k+* has to be sufficiently smaller than k* such that condition 
number 2 as set out for eq. ( |T6"D is satisfied. 



The constants C\ and C 2 should be 



k+ _, rlh 

A'pp(k)k~ 1 — 



Co 



k* ~ rib 

A 'pp{k)k- Z — 



Assuming some values for G\ y C 2 and the starting wavenumber k±+, eq. (|T7]) can be used 
to obtain A'P p {k = ki) for successively smaller fcj's (at, say, evenly spaced intervals). The 
method adopted by Croft et al. ( |1998| ) is equivalent to keeping only the first term within 
the square brackets on the right hand side of eq. (|i7|) i.e. a simple differentiation. 

Because we are interested only in small fcj's, and because P p (k) generally falls rapidly 
with increasing k, especicially at high fc's, we can choose k+* to be some sufficiently large value 
and simply set C\ = C 2 = 0. We will see that our method is robust enough to consistently 
yield good agreement with the input power spectrum on large scales (small ki), even though 
one is making an error on small scales by approximating C\ and C 2 as zero. Also, strictly 
speaking, k++ should be chosen to be smaller than k+, which is the k-value above which the 
distortion kernel is no longer described by the linear kernel in eq. Jl5] (see also conditions 
for eq. [0). We will not be careful about it, and will see that one still obtains the correct 
P p on large scales, again because P p falls rapidly with k. 

For clarity let us call the method of simple differentiation following Croft et al. ( |1998 ) 
method I, and the alternative that we propose here method II. 

Fig. |l| shows a comparison of the inversion using method I versus II. Since it is generally 
difficult to judge differences between power spectral shapes in a log- log plot of the power 
spectrum, we instead show the fractional error in the inverted power spectrum as a function 
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of k. The input mass power spectrum is that of a SCDM (Standard Cold-Dark-Matter) 
universe, with Q m = 1 and h = 0.5. We have chosen the parameters 7 = 1.5, fcjj = 
0.11(km/s) _1 and kp = 0.12(km/s) _1 in the input W^ p (eq. [0]). The latter two values 



correspond to the choice To = 10 4 K (eq. JTT] & [T^]). (We will later alter k^ to mimic 
the effect of nonlinear redshift distortions.) The length scales shown correspond to those 
considered by Croft et al. (|1998|) . Since we are only interested in shapes here, the power 



spectra are normalized to agree at k — 0.005(km/s)~ 1 . It can be seen that method I gives 
an inverted mass power spectrum which, on large scales, is systematically less steep, or less 
steeply falling, than the input (i.e. the slope of the inverted power spectrum is less negative 
than the actual one; see Fig. [3| for a log- log plot of the power spectra). 



To understand this result, let us go back to eq. ( |Ty|) and rewrite it as follows, approxi- 
mating C\ and C2 as zero and k^ as effectively infinite: 



2tt dP f 



= A'P p (h) + E{h) (19) 

h' 1 



(l + !3f) 2 k, dk, 



On the left hand side of eq. ( |19|) is essentially the estimator of Croft et al. (|1998| ; method 
I) for the shape of the mass power spectrum. The first term on the right is the true mass 
power spectrum, and E is the error of method I. Approximating a realistic power spectrum 
P p (k) as k~ n [l — t{k)] where e(k) is positive and is an increasing function of k, it is not hard 
to show that d(E/ P p )/dki > by expanding to first order in e, which means E is decreasing 



with ki slower than, or at most as fast as, P p is. Eq. (|T9"D then tells us method I would 



systematically give a flatter estimate of the mass power spectrum than the true one. The 
limiting case of d(E/P p )/dki = occurs when the input power spectrum P p obeys a strict 
power-law, in which case E has the same shape as the input, and method I recovers the 
shape of the true power spectrum. 

It can also be seen that both method I and II fail on small scales (large k). This should 
come as no surprise because no attempts have been made to model the small scale effects 
in the inversion procedure laid out in eq. flT7D (see also eq. flTE|). We "pretend" that we do 
not know the actual full distortion kernel (eq. [fRJ]), but instead assume only knowledge of 
the large scale distortion kernel (eq. fll5"| ) when carrying out the inversion. Moreover, for 
method II, we have not been very careful in selecting the value of the constants C\ and C2: 
we simply set them to zero. 

To estimate the effect of nonlinear distortions on our inversion procedure, we decrease 
kl to 0.028(km/s) _1 , and show the outputs of method I and II in Fig. 0. Here we are taking 
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advantage of the fact that the factor of exp[— {k\\/k^) 2 ] in eq. ([14]) is commonly used to 
model nonlinear distortions in the case of galaxy surveys (see e.g. [Peacock fc Dodds 1994| ). 
By raising the scale of nonlinear distortion by a factor of about 4, we hope to gain an idea 
of how the as yet poorly understood nonlinear distortions on small scales might affect the 
inversion of the power spectrum on large scales. However, it remains to be checked using 
simulations how realistic this choice of scale, or this particular parametrization of nonlinear 
distortions, is. The agreement on large scales for method II is not as good as before, but 
is still within about 7%, and is better than that of method I. It is possible to improve the 
agreement by playing with the input parameters C\, C2 or k^. We will not pursue that here. 

We also show in Fig. [3] what the inverted and input power spectra look like in a log-log 
plot. The subtle differences in the shapes of the power spectra are still possible, but harder, 
to discern in such a plot. 

Lastly, we have assumed in all tests above that the input (3f (eq. fl5fl ) is known when 
performing the inversion. In practice, there is an uncertainty due to the lack of knowledge 
of the precise values of 7 and In Fig. |], we show a case where the inversion 7 is chosen 
to be slightly different from the known input 7. (The actual value for 7 in the real universe 
is likely to have a narrow range 1.3^7^ 1.6 c.f. |Hui fc Gncdin 1997] ). The impact on the 
recovery of the power spectrum shape appears minimal, for the small change in 7. Similary, 
at a redshift of around 3 should fall in a narrow range (close to 1), for reasonable values 
of Q m and the cosmological constant today. 



5. Conclusion 

The main aim of this paper is to draw attention to the fact that redshift distortions 
in general make the three-dimensional redshift-space mass power spectrum anisotropic, and 
so the inversion from the projected one- dimensional redshift-space transmission power spec- 
trum to the three-dimensional real-space mass power spectrum involves more than a simple 
differentiation (eq. [0]). Given a kernel (W^ p in eq. ||) that relates the relevant power 
spectra, it is possible to perform the inversion by essentially inverting a triangular matrix 
proportional to W fp (eq. & [§§). 

We have demonstrated this idea with a kernel W^ p that is motivated by linear pertur- 
bation theory. In general, we find that a simple differentiation method tends to make the 
inverted three-dimensional real-space power spectrum flatter than it really is on large scales. 
A procedure that remedies this is outlined in eq. fll7|) . We have referred to the former 
(straightforward differentiation) as method I, and the latter as method II. 
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As we have remarked before in §C2, if the three-dimensional real-space mass power 
spectrum P p (k) obeys a strict power-law, the simple differentiation procedure of method 
I will recover the same power-law. That we observe a deviation of the method-I-inverted 
power spectrum from the input in Fig. |T] to |] is a reflection of the fact that a realistic power 
spectrum (such as CDM) is only well-approximated by a power-law for narrow ranges of k. 
The differences between method I and II are not easy to discern in a log-log plot such as the 
one shown in Fig. especially when one is dealing with observed or simulated data with 
noise. The systematic error of method I is nonethless present, and more easily seen in plots 
such as the one in Fig. |2|. A cursory inspection of some of the log-log plots of the inverted 
power spectrum in Croft et al. ( |1998| ) seems to indicate that method I does give a slightly 
flatter power spectrum compared to the input on large scales, but it certainly should be 
more carefully quantified with simulations. 

An interesting consequence of redshift-space distortions is that the one-dimensional 
transmission power spectrum P* is no longer guaranteed to be monotonically decreasing 
( [Kaiser fc Peacock 1991| ), unlike in the case where the distortion kernel is trivial (i.e. set 



W$ p to constant in eq. |13|). Method I, where P* is simply differentiated to obtain the mass 
power spectrum P p , could then give a negative mass power spectrum and fail dramatically 
(see eq. [[19[], and discussions that follow). 

For a power-law mass spectrum with P p (k) oc k~ n , it is simple to show using the linear 
distortion kernel (eq. pj), together with eq. ([13]), that dP f /dk {l > if [-J]P f + 6/3/ + 1 - 

< n < [y/pj + 6/3 f + l-(l-^/)]/(l+/3/). For the parameters Q m = 1 and 
7 = 1.6, which were adopted by Croft et al. ( |1998| ), this implies a range of —1.62 < n < 1.17, 
within which method I should give a negative inverted mass power spectrum. Of course, a 
strict power-law power spectrum with n in such a range is not very interesting because 
it gives a diverging one-dimensional transmission power spectrum ( n > 2 is required for 
convergence). It is for this reason that a direct comparison cannot be made between the 



above simple estimates and a test-case shown by Croft et al. ( 1998 ) in which the initial 
power spectrum has n = 1. In fact, for a scale-free initial power spectrum of k~ n , the stable- 
clustering mass power spectrum at high k should be & _6 /( 5 - n ) (|Peebles 1980|). For n — 1, 



this means the high k mass power spectrum asymptotes to k 1,5 , which is still not enough 
to regularize the integral for the transmission power spectrum (eq. |13| together with eq. 



|15|). It is likely that in their test-case, an effective ultraviolet cut-off is imposed by the finite 
resolution of the simulation, or else the true nonlinear redshift-distortion kernel provides an 
effective regularization. Note, however, that there is no small-scale Jeans smoothing to help 
here because Croft et al. ( |1998|) used N-body-only simulations. 

Nonetheless, we should emphasize that while, on large scales, the direction of the sys- 
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tematic error of method I is on firm ground (i.e. systematic flattening), the size of the 
error, on the other hand, is subject to further investigation. We have used linear theory 
in §^ to estimate the magnitude of this effect. The linear calculation gives us the effective 
"bias-factor" (let us call it bp, which equals 2 — 0.7(7 — 1) in eq. |T5|| ) associated with the dis- 
tortion parameter /?/ (i.e. /3f = fa/bp), but this might not be the correct bp in the presence 
of nonlinear fluctuations on small scales. 

For instance, according to arguments by Scherrer & Weinberg (|1997|) , the galaxy-bias 
b g (different from bp above, see below), defined as the ratio of the galaxy to mass two-point 
function on large scales (where galaxy is related to mass by a local transformation), generally 
involves all the higher derivatives of the local transformation around 5 = i.e. not just the 
first derivative, as is in the case of a linear perturbative calculation. 

One might naively think that, in the case of the Lyman- alpha forest, bp should be set 
equal to the analog of the galaxy-bias defined above, namely the ratio of the transmission 
to mass power spectrum on large scales. Let us denote this latter quantity also by b g . This 
number is smaller than 1 in our case (because the exponential (e~ T ) suppresses large fluctu- 
ations), which means the redshift distortions are more pronounced than we have assumed in 
§^j, and the systematic error of method I should be even larger! 

However, there are at least two reasons to suspect that this is not the correct conclusion. 
First, as emphasized recently by Dekel & Lahav ( |1998| ), the bias factor that shows up in the 
linear distortion kernel (bp) is not necessarily the same as b g defined above, because of the 
nonlinearity of the biasing transformation (see below for a discussion of the transformation 
relevant in our case; we do not suffer from stochastic biasing, however). In fact, according 
to Dekel & Lahav, even the form of the linear distortion kernel could be slightly modified 
(see also [Pen 1997| ). 



A second, perhaps more important, reason is that the mapping from the mass density 
field p to the transmission / actually involves two local "biasing" transformations with the 
redshift-space distortion in between. First|], the mass density in real-space is related to the 
neutral hydrogen density in real-space through the local transformation n m oc p2-o.7( 7 -i). sec- 
ond, this is then "distorted" into the optical depth in redshift-space through r oc [p 2 ~°' 7< ' 7 ~ 1 * ) ], 2 , 
where [] z denotes a quantity in redshift-space; finally, another local transformation maps the 
optical depth in redshift-space to the transmission in redshift-space: / = e~ T . It is plausible 
that the first two steps determine the correct value of bp, while the last merely shifts the 
overall normalization of the final redshift-space power spectrum of the transmission on large 



scales (see [Hui 1998|) . 



D We ignore the spatial dependence of the thermal profile to simplify the discussion here; see eq. 
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Lastly, there is of course the possibility that even the form of the large-scale distortion 
kernel in eq. ( |i~5|) might be incorrect, or as is more likely the case, that the linear distortion 
kernel applies only at (a small range of) the largest scales (see e.g. |Uole et al. 1994] ). 
This is related to the question of how nonlinear clustering on small scales, or translinear 
clustering on intermediate scales, affects the large scale behaviour of redshift distortions. 
Couple this with the effects of nonlinear local transformations: we indeed have a complicated 
problem here. Some of these issues are begining to be addressed (see [Fisher fc Nusser 1996| ; 



[Taylor fc Hamilton 1996| |Scoccimarro 19981 ; |Hui et al. 1998| for the former, and (5cherrer~S; 



Weinberg 1997| ; |Dekel fc Lahav 1998j ; [Hm 1998| for the latter). The good news, at least in the 



case of the Lyman- alpha forest, is that the form of the relevant local "biasing" transformation 
is known exactly. 

Both issues, the problem of translinear or nonlinear redshift-space distortions and the 
problem of biasing in real as well as in redshift-space, are obviously of great interest in 
the wider context of large scale structure and galaxy surveys. Analytical calculations, with 
reality checks using simulations, would be necessary to address these questions. We hope to 
pursue aspects of these issues in a future publication. 

As this paper was nearing completion, the author became aware of a preprint by Nusser 



& Haehnelt ( [1998|) who also considered the effects of redshift-space distortions on the for- 
est, not in the case of inversion from the transmission power spectrum to the mass power 
spectrum, but in the case of the recovery of the mass density field itself. The author thanks 
Rupert Croft and Albert Stebbins for useful discussions. This work was supported by the 
DOE and the NASA grant NAG 5-7092 at Fermilab. 
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Fig. 1. — A comparison of the Croft et al. (1998) inversion (method I; dotted line) versus 
inversion set forth in eq. flTTD (method II; dashed line). On the y-axis is 5P P /P P where 
fipp = P p utput — P p , and P p is the input real-space mass power spectrum, and -Poutput is the 
output. The input power spectrum is that of the SCDM model, and the input parameters 
are 7 = 1.5, k?< = 0.11(km/s) _1 and hp = 0.12(km/s) _1 (eq. ||14|| ). The inversion parameters 
are k M = 0.14(km/s) _1 , G\ = C 2 = (eq. ]I7J). The inversion j3f is assumed to be the same 
as that in the input (i.e. 1/1.65, see eq. KJ). The outputs are normalized to match the 
input at k = 0.005(km/s) _1 . 
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Fig. 2.— Same as Fig. | except that fcjj = 0.028(km/s) _1 . 
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Fig. 3. — The solid line is the input linear power spectrum. The dashed line is the output by 
using method II (eq. 1171), and the dotted line is the output by using method I i.e. simple 
differentiation. The parameters are the same as in Fig. [| 
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Fig. 4. — Same as Fig. |T| except that the inversion 7 — 1 is chosen to be 0.3 instead of the 
input value 0.5. 



